A Relativistic Flux-tube Model for Hybrid Mesons 
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A number of authors have considered potential models for hybrid mesons. These frequently involve 
approximating the vibrating flux-tube by a set of beads, and making an adiabatic approximation 
which gives rise to a static inter-quark potential which has an effective repulsive - potential. We 
show that this approximation is almost certainly wrong. Since the beads are presumably massless, 
the correct approximation requires the solution of a Klein-Gordon-like equation treating the beads 
and the quarks on the same footing. We show how to solve this in the one-bead case for massless 
quarks, and And the spectrum has an unexpected degeneracy. We generalise this to the As-bead 
case, which can still be solved exactly for massless quarks, and show how to renormalize the energy 
to obtain a plausible spectrum. We give a generic method to solve the equation for massive quarks, 
and use this to derive a different non-relativistic equation. The results show a smooth behaviour 
both with respect to the quark mass and the number of beads. 

PACS numbers: 



I. INTRODUCTION 

Many authors have considered models for hybrid 
mesons. In most cases, the method used is similar in 
spirit to the original work of Isgur and Paton . In its 
simplest form, the meson is regarded as a qq state with 
an interaction potential produced by an oscillating string 
or flux-tube. A physical model for this consists of mas- 
sive beads linked by a linear potential (Fig [TJ) These are 




FIG. 1: A hybrid meson model 

treated non-relativistically in an adiabatic fashion and 
the oscillations are restricted to being perpendicular to 
the line joining the quarks. This basic idea has provided 
a rich source of models. There are justifications for var- 
ious models from lattice-gauge theory (0,0,0). More 
sophisticated analyses include non-adiabatic corrections 
0. and a number of phenomenological analyses: see 
s~K-(0! 0)- The work by Barnes, Close and Swanson 
[8( discusses a non-adiabatic model which has a similar 
starting point to the discussion below. 

This paper discusses a formulation which escapes from 
some of these restrictions. Firstly we show how to elim- 
inate the centre of mass from a N-bead non-relativistic 
model. We use this to examine the adiabatic approxi- 
mation and show that it can lead to highly misleading 
results in a very simple case. Then we write down a 
model where both the quarks and beads are treated as 
Klein-Gordon particles. The resulting equation can be 



exactly solved in the case where both quarks and beads 
are massless, and we write down a variational expression 
for the case where the quarks are massive. There are 
a number of subtleties in this model and the spectrum 
has an unexpected form. We have not made an attempt 
to turn this into a realistic model for conventional and 
hybrid mesons. 



II. CENTRE OF MASS DIAGONALIZATION 

The assumption that the quark motions are indepen- 
dent of the flux tube is not particularly logical for small 
quark mass. We therefore consider the extension of the 
flux-tube hybrid model to include the quarks. If we as- 
sume the quarks and the beads have small masses, we 
can write 
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, fi N = m N 
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where the flux tube mass m / need not be a constant and 
we are allowing for the possibility of the quarks having 
different masses. 

The full Schrodinger equations is 



1 



2mi 



■V? 



1 



JV-1 



2mN 



-V 



N 



k=2 



o-l N (r, y n ) * = E^> 



In the spirit of the adiabatic approximation, we could 
assume that the beads are restricted to oscillations per- 
pendicular to the tube, so V| is 2-D, but and Vjy are 
3 dimensional, so the quarks can move in 3-D. 

If the quarks are light, they need to be treated on the 
same footing as the beads, which means we must separate 
the centre-of-mass motion from the relative motion terms 
in the sum over the KE. The CoM given by 
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Writing 



N 



h = R = 



-n-i(i>2) (1) 



eliminates any cross terms of the form. ^ J^- Note that 
this reduces to the usual expression for the 2-body case. 
This can be inverted, via 
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n-i =n-Ti 

This gives 

N 1 i N i 

n— 1 n— 1 

We now want to diagonalize the KE operator. Write 

n— 1 

where the matrix G is symmetric and tri-diagonal. 
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We can then diagonalize this using 
u = Sf, S T GS = A 
which allows the kinetic energy to be written 

i d 2 v 



2mn dr 2 



T— v 2 

n— 1 



It is then possible to solve the equation analytically if 
we make two further assumptions: that the motion is 
adiabatic, so that the beads oscillate perpendicular to 
the line joining the quarks and that P.E. 

N N r o-il/2 



Hi 



(*)' 



In {r,y n ) =o-^|f,-| = <t^2 

i=l i=l 

(where yi is perpendicular to the line joining the quarks) 
can be expanded: i.e. 

N 



V N {r,f n ) + — 



2r 



so the flux tube interaction gives rise to an oscillator, but 
with a strength that increases with N. 

N 

Since u = Sf , the P.E. term can be written [ff 1 = 

i=l 

N 

\uf] This allows the N-bead equations to be separated 

i=i 

into N uncoupled equations, where the effective spring 
constant is given by 

K eff ~ - 



III. FAILURE OF THE ADIABATIC METHOD: 
ZERO-BEAD CASE 

We would like to check the validity of this expansion in 
a very simple example: the bead case which be solved 
exactly in both 3-D and 2+1 D methods. The basic equa- 
tion is 



(3) 



2m R 



V'-ar + E nl vj nl (r) = 



1. Exact solution 

As is well known, for 1 = 0, after conversion to a di- 
mensionless form via 

nt \ V n a f2ma\ 1/3 

y = p{z - z n ) ,z n = — ,p= l-^ - ) 



one obtains the Airy equation giving 

E n - 



2m 



Vn 



where y n is the n-th zero of the Airy function 



2. Conventional Variational Solution 

We can, of course, solve ((U by a standard variational 
approach with a Gaussian- wave function. The variational 
energy 

h 2 „ 3 2a 



E=—0 2 - 



gives 



2m 2 /Vtt 



/? = 0.867 (^) 1/3 ,i^ 1.865 AV2X V3 



3. 2+1-D solution: 

We can now solve the 3-D equation by a "phonon" 
approx. This means put r — (x,y) so 

\r\ = \J x 2 + y 2 « x + — 
and the trial wave function is 

0„ (x,r) = e-^"*) 2 e~^^T 
The 2-D "phonon" equation has the form 

/ fi? ( d 2 | 1 d £ 2 \ ay 2 | E \ Q 
\2mn. \dy 2 y dr y 2 J 2x ™ / n 
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with the ground state solution 

„, 1 lam „, , , ^ / a 
n V x V mx 

This then gives the corresponding 1-D equation 

\ Zm ux V mx J 

with a variational solution 

2m 2 V^/ 3 V^r 
which can be minimised numerically for /3 and E. 

4- Simultaneous minimisation 

To justify a technique we will use later, we can min- 
imise over both /3, f3' simultaneously. This means doing 
a non-linear minimisation for both the phonon and lon- 
gitudinal modes: i.e. split the equation into transverse 
(phonon) and longitudinal motion as before 

(-& + ^ + + A> (x,y) 

= E<j) (x, y) 

and using a trial wave-function 0o (^j y) = 
e —{0x) /2g— (@ v) /2 ^ perform simultaneous non-linear 
minimisation on /3, /3' 

5. Results 

The results, using m g = 300Af eV, <r = ISOOOOMey 2 , 
are shown in table HI We conclude that the phonon ap- 

TABLE I: Comparison of zero-bead results 

ft Energy 

3-D exact 1.3203 1113.45 

3-D variational 1.3850 1116.63 

2+1-D phonon approximation 0.89668 1567.87 
2+ 1-D minimization 1.3846 1116.62 

proximation is disastrously bad. It can be argued that 
the model is unrealistic, since the transverse motion could 
not be regarded as adiabatic in this case. However, we 
find similar results for a one-bead example, but there are 
no exact results to compare it. 

IV. KLEIN-GORDON EQUATION 

The assumption that the quarks can be treated non- 
relativistically is not sensible when the quark mass is 



small. Further, it is hard to justify using a non- 
relativistic approximation in the limit that N — > oo since 
the bead mass goes to zero and the Schrodinger equation 
becomes meaningless. Hence we would like to have an 
equation which works for both massless beads and mass- 
less quarks. 

We consider a Klein-Gordon equation for the both 

{pic 2 + mjc A ) i/f = 0, mi - m + -± (4) 

where we have assumed that the confining potential 
transforms as a Lorentz scalar (assuming it transforms as 
the 4th component of a 4- vector leads to other problems: 
see H). Obviously a realistic model would require a Dirac 
equation for the quarks. 

The spectrum should be well-behaved and independent 
of the mass in the [1q — » limit, since the K-G equation 
becomes 

(-^ 2y2 ~ ^ + V2 ) ^ = ( 5 ) 

We need to reduce the n-body K-G equation to a 
Schrodinger-like equation with no relative time. This 
does not appear to be in the literature. One cannot use 
a conventional reduced mass reduction: see ([HI)- 

If we have n particles with individual momenta pi , we 
can obviously write the total momentum 

which is correct relativistically or non-relativistically. 
Each bead or quark satisfies 

(p 2 c 2 + m|c 4 ) ijj = e 2 -0, m, = M + 

so by analogy with the preceding we write 

Pi = + qi - q%-i 

Here and in what follows we use E to refer to total system 
energies and e to refer to individual component (quark, 
bead or phonon) energies. This is specified more carefully 
in Appendix A. 

Then the <?'s can be reconstructed 

qi=Pi-%P, 

k 

k £ a 
qk =Pk + qk-l - %P = Y,Pi e~ P 

Note that qk is still purely space-like: and the total 

n 

energy E — ^ Si . Effectively in this process we are 
i=i 

replacing the masses of the non-relativistic particles with 
the corresponding energies. 
Hence 
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giving 



p2 

E 



E t 



(6) 



-h 2 c 2 



E 



(V ?; 



+ V 2 )xJj = Mc 2 tp 



so the overall K-G equation becomes 
£i(p 2 c 2 + m?c 4 )* 



E 



(7) 



This equation is correct but apparently useless, since 
the individual £i are unknown: in the equivalent NR 
treatment, the energies are replaced by masses whaich 
are, of course, known a priori. Below we show how to 
find the Si by a variational trick. 



V. MASSLESS QUARKS 



The K.E. can be diagonalized with eigenvalues 
l/e+, 1/e_ given by 



E = 2ei 



■£2, £+ = £!,£- = - ^ 



_1_ 

E2 



2E1+E2 



w± = 7f (ft ± h) ,r 2 +f 2 =u 2 + 



2f 2 f 1 



(v 2 



(note that this is just the relativistic analog of the re- 
sult in Merlin and Paton [f|. Hence the K-G equation 
becomes 







° 2 r 






+ -. 
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£_ 



1 



■0 = £0 



We consider in turn the (trivial) zero-bead solution, 
the one bead solution and show how the solution can 
be generalised to the iV^-bead solution. Then we show 
how the -/Vb-bead solution can be renormalized to give a 
sensible spectrum. The starting assumption is that there 
is a scalar confining potential, V = or equally divided 
between the quarks and the beads, so that (U|) becomes 



™2 4 2 4, 2 , 



A. Massless zero-bead solution 

The energy equations simplify to 

, Mc 2 M 

E = Mc 2 ,E R 



Mc 2 



Then K-G eqn becomes 



which separates into two phonon equations: 



2^,2 



4^ 



-0 = e±ip 



Note the distinction between the e, (individual energies 
of the quarks and bead) and e± (individual energies of 
the phonon states). This has obvious solutions: 



$± - 2S 

Eq = i+ 



1/2 



'•/?•) 



2v^+e 

- i 3cha 

■ S— — £-l -f" fc_ — 



= (2n± + L + |) 



To solve this, we will assume that the quarks carry frac- 
tions cti of the total energy of the system, and regard a.\ 
as a variational parameter, so 

ei = a.\E = e 3 = £+, £2 = ob£, «b = (1 - 2ai) 
£_ = £ 2 (ai) - , 3cto 



E -V 2 + -^)r Er 



with an immediate solution 



/? 2 = — ,E 2 = Qhca 
2nc 



gives a minimum value of a\ = 1/3: i.e. the quarks and 
the bead carry an equal fraction of the energy, which is 
what we expect. This gives 



2ch 



1/2 



2ch 



(3) 1/2 ,E 2 



27 cha 



B. Massless One-Bead 

The interparticle potential becomes (f„) 
§ (ft + ft-i)- Hence the sum over potentials is 



V 2 = 



x (- + — 



E3 

2f Q f 



Physically we can interpret the two solutions as show in 
fig[U the + solution corresponds to the two quarks mov- 
ing in opposite direction while the bead remains station- 
ary at the centre of mass, while the - solution corresponds 
to the quarks moving in the same direction but opposite 
to the bead. 

It is worth noting the peculiar feature of the solution: 
the two phonon states have the same excitation spectrum 
but different (3's. 
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FIG. 2: Two modes for one-bead model 

C. Massless n-bead solution 

The overall K-G equation j7]) becomes 



E 



2 „2 



(gj — g£- l) C 



Ei 



If we assume that the e^s are equal 

E 



N + 1 



(8) 



then this can be solved (we justify the assumption below). 
The previous method generalises: 



h c N 



f^-E(n + ri_ 1 ) 2 W' = 



£ 4 



where 



G 



2 -1 
-1 2 -1 

-1 2 -1 
-1 2 



2 1 
1 2 1 
1 2 1 
1 2 



and 



= Sf,S J GS = A 



(9) 



The structure is obviously similar to 

The KE term can be diagonalized as before. The same 
operation that digaonalizes G also diagonalizes G' . The 
eigenvalues are given by [III, [HI- 



A' l N = [ 2 sin 



2(JV + 1) 



This then separates to give a set of N uncoupled 
"phonon" equations of the form 

/ h 2 c 2 N Na 2 Ajv , 2 \ , . , 

( — Af V 2 Ui + ——A%_ i u 2 J i> = W 



which have the usual SHO form, with different interpre- 
tations of the energy and potential. With a modification 
of the usual substitutions 



Ajv_ 



cat ( m 1 _2_ 

Ai 2h ~~ CULj \ N+1 ) 2h 



(m 2 - 

= (2n 2 + k + §) sin (^) 



(10) 



we can write the solution as 

The energies can be summed to give the total energy 



N 



i=l 

3 chajN+V) I-'" 



( Jv + i ) 



N 

pTV _ -AT _ 3 c7to(iV+l) ^ „;„ / 7ri 
— fc i — 2 2E« Mn \ 7V+1 

a=l 



2E- 



l-oos( 1 ^ T ) 



so 



There are several problems with this solution. Not 
surprisingly, the total energy is divergent as N — > oo. 
Secondly the curious degeneracy noted above remains: 
the energy of the n th and (N — n) th levels are equal, even 
though they are described by different values of j3. This 
violates one's intuitive idea that large corresponds to 
small E and vice-versa. Finally, the individual energies 
of the phonon modes are 



3cha tt 1 / 2 



N+l 



3cha 7r 3 / 2 i 
2 2(AT+1) 



so the excitation energies are vanishingly small in the 
N — > oo limit. 



D. Energy Renormalization 

In order to obtain finite answers, we subtract a energy 
AE N from the total non-renormalized energy E N given 
by (fTT)|) to give a (fixed) renormalized energy E° for the 
ground state: 

E° = E N - AE N 

By writing 

N 

AE N = Y,if (0,0)- J° 

i=l 

gives the renormalized energy of each phonon mode as 



cha (N + 1) 
2E° 



{2rii + li) sin 



N+l 



(11) 
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This now has a sensible limit for all modes in the N — > 
oo limit 

eTK,/i) = ^ ) {2n l + l l )m 

It makes sense to have the energy fixed by the zero-bead 
solution, so 



The second term, in particular, causes considerable trou- 
ble since it mixes the phonon states. We provide a very 
rapid numerical algorithm for taking this into account. 
In principle one can include a bead-mass: there is no 
need to do so, and we show that one can in fact find 
a non-relativistic reduction of the equation even with a 
massless bead. Again, we analyse the 0-bead, 1-bead and 
n-bead solutions. 



E° = VGhca 



(12) 



and the total energy of the excited phonon states is 
then given (in the N — > oo limit) by 



E°° (n k ,l k ) = V&dm (l + ^ ( 2n * + *0 *J 



(13) 



VI. MASSIVE QUARKS 

We have ignored the quark masses: this introduces two 
extra types of term, via 



2 4 2 4 i 2 i 

m i c — [i i c + ^crrc + 



A. Massive quarks, zero-bead case 



We consider the general unequal-mass case. The K-G 
equation is 



h 2 c 2 If 1 
V 2 + — ^?c 4 + fi lC 2 ar + -a 2 r 2 



j- Ui 2 2 c A + [i 2 c 2 ar + iaV j ) r = Ei 



In order to solve all of this variationally (and to connect 
this easy example with later equations), we proceed as 
follows. 

1. Assume that the quarks carry fractions a\, a 2 of 
the total energy of the system. Then 

E\ = a\E, £ 2 = a 2 E, er — a\a 2 E, u\ + a 2 = 1 

obviously < a±, a 2 < 1 

2. Hence the K-G equation becomes 



Non-relativistic Reduction 



(_ RV V 2 + + + c 2 0~T ( i± + 

I QlQ2 Ql Q2 y«l 

= E 2 tp 



4> 



(14) 



3. Use a Gaussian trial wave- function, as suggested by 
the massless case, giving the variational energy: 



y> _ h 2 c 2 3 a2 , fJ-lc 4 

Ot\Ot2 ^ rVl 



E 2 = 



c 2 a (m. + uz) _L_ + si 

L " \ ai ~ a 2 J fj^/jf T 4 2ai 

This obviously contains the massless case. 



3 



(15) 



It is not obvious that this technique gives a sensible 
non-relativistic reduction. If the masses are large and 
equal, M = 2/i + ^g- where Eb is the binding energy. 
Putting 



E = Mc 2 = me 2 + fi 2 c 2 + E 



£i fts nic 2 ,s 2 m M2C 2 ,£fl 



B, 
M1+M2 



Ignoring terms with V 2 , .Efj in (fl"4"l) gives 



V 2 + err - E B W> = 



which is the Schrodinger equation. 

However, we wish to use a variational method, and 
it is not conventional to find the reduced mass by this 
technique! We get a mass term from (|T5|) of the form 

( 2 4 2 4 \ 
Mi£_ _|_ £a£ £; 2 ) . In the infinite mass limit, this 
CKl CK2 J ' 

gives an apparent mass term 

m 2 = (/i 2 (1 - ai) 2 + [i\a\ - 2/n (1 - ot\) /iaail 
= (/ii (1 - ai) - Ai 2 ai) 2 
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Minimising this gives (obviously) a-y 
duced energy" becomes 



mi 



so the "re- 



£r = a\aiE = urc 



as expected. 



2. Coulomb-interaction and non-zero angular momentum 



This introduces two extra terms. There is a constant 
term with same 



d = 



n 2 r 4 



as in (flUl). The second term is 



(17) 



We can generalise this by including a Coulomb-type 
interaction (see Olsson [13]) V{r) = ~ which is the fourth 
component of a 4-vector, so that 



£l -Yr) =£l \ £l --r 



ei4r 2 



This adds a term 



2nE 



to the K-G equation. 

We can also generalise this to any angular momentum 
I: both these generalisations require us to use a modified 
trial wave function of the form 

H 2 r 2 

We will write the variational energy in the generic form 

E 2 (m, n 2 , £) = h 2 c 2 d^ 2 (aO P 2 + nEd-x {(3) + 
c^do ( Pj , ai ) + oc 2 d x (/?) + 

where 



d-2 
d-i 
d = 
di = 



90 r( P +i) 



4i(i+l)+4p+3 



Hi. 



do = 



M2 
a 2 



(16) 



r( P +2) 



This generic form is useful since it generalises to the 
one-bead and n-bead case, and it can be solved by a 
rapid algorithm which is outlined in appendix B. For the 
purposes of renormalizing the n-bead energy later, we 
regard the zero-bead energy as the renormalized energy 
of the ground state 



E u 



in the sense of (CE 



E(ni,^2j) 



B. Massive quarks: one-bead 

To solve the unequal-mass one-bead case, we will as- 
sume that the quarks carry fractions a\, 013 of the total 
energy of the system. 

£1 = aiE, e 2 = a B E, e 3 = a 3 E 7 as = (1 — a% — 0-3) 



2c z 



si 



£3 



a 3 



E V ai 



(18) 



is more difficult to handle 

In the equal mass case, (fTTf becomes: 



do = 2^ 
a 



and the linear term (fT8l) becomes 



V2 



ficr 
a\E 



+ t t 2 _| 1/2 (l + 0(cos 2 (9)) + ...) 



Ignoring the higher order terms gives 



h 2 c 2 ( v 



2fi 2 c 



2„4 



a ^ E \ V2c 2 na \u 2 



-2 I 1 / 2 



r -2 
+ 



tp = EtJj 



The linear term (|18p gives us a new term in the varia- 
tional energy 



A 2 V2c 2 ^ J J e - +<e-^- u - \u\ + u 2 _\ 1/2 ul uidn _J« 



V2c 2 ^tio(P+,P-) 



which can be written as d\ (a, (3) — V2c 2 ^fo (/?+, (3-) 
in the notation of (fl6|) . We show how to evaluate /q 
in Appendix C. . Hence the overall variational solution 
becomes 



£2 = 3^ (/3 2 



Hi 



2^c 



2^,4 



^ 2 ^/oo(/3 + ./3-) + §£(^ 



/3i 1 /3? 



which can be rapidly solved by the technique of Appendix 
A. 



1. Non-relativistic reduction 

It is useful to have a non-relativistic reduction of this 
equation. The conventional reduction of allowing the 
bead and quark masses to become large will not work, 
since we believe the bead mass is zero.. We need to re- 
duce (|VI Bj) to a corresponding NR equation for compar- 
ison. Putting E = 2/ic 2 + E B gives 

a x E 2 - 2 fJ , 2 c i w 2/i 2 c 4 (2ai - 1) + 4q!i^c 2 .Eb 
= Active 2 Eb — 2/i 2 c 4 as 
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The K-G equation then becomes 



Finally the mass term gives 



= E B i> 



vl 



1 9 1 a 1 ) a O 1 1/2 



so the bead has acquired a pseudo-mass aiia^A 4 
This has the variational solution 

F 3 h 2 f 0-\ I ff— \ I 

^7-/oo (/3 + ,/3-) + ^c 2 |f 

Note that the masslessness of the bead does not give rise 
to a Schrodinger equation with zero mass. The usual 
limit of a = 1/2 which would minimise the mass term 
does not minimise the energy. 

C. Massive Quarks: n-bead case 

We want N B = N — 2 beads and 2 quarks, with dif- 
ferent masses, so that we have N particles in all and the 
quarks have a fraction ct\ , on of the total energy, so each 
bead has a fraction 



OLB 



1 — OL\ — Ct2 



(19) 



N -2 

of the total. Now the K-G eqn becomes 

v -h 2 (v, - v.-i) 2 c 2 + + v % ) 2 

The KE and V 2 terms give expressions similar to ^ with 
G and G' are now given by 



G 



G' 



a\ + 1-1 



-1 2 -1 
-1 2 



-1 

-1 a' 2 + 1 



1 1 



2 1 
1 2 1 



1 a' 9 + l 



with 



1 — Oil — OiN Ob 
Oil 



1 a x {N-2) 
The term in fiV becomes 
2 / HxVi (J-nVn 



2c' 



El 



ocr_ I Mi Fi| 
E 



oi\ 



Fat I 
a N 



(f 
















V ai 





as before. There are several new difficulties: 

1. for non-equal masses, diagonalizing the KE matrix 
does not diagonalize the PE matrix. 

2. the linear terms like /ii |ri| involve all the Uj's, so 
the fiV is very complex. 

3. there are (N^ + l) different (3's (one for each phonon 
mode) which must be found variationally. 

4. in addition, we have a\ and a^v which must be 
regarded as variational parameters. 

Surprisingly, this can be approximately solved by a quite 
simple numerical technique: the algorithm is as follows. 

1. Guess a\,ajsi 

2. Diagonalize G and find eigenvectors e\ 

3. Ignore off-diagonal terms in e T G'e = A' 

4. Write T\ — Y^ e j n j so that 



ri 



,(i)Ji), 



ej 'e' k 'UjUk 



5. Then (approximately) the expectation values re- 
quired for the linear interaction can be found (See 
Appendix D): 



/i(/3i) = <N> 




1/2 



6. This gives us a total (non-renormalized) energy of 
the form: 

^=fiVEA^(l j + |) + ^ + f 

w {^h (A) + m) + ^e|o, + |) 

to be minimised to find the /3j. Obviously if \x\ = 
/ijv = this separates: it still has the form of IB11 



7. After finding the energy f3j. the energy must still 
be minimised by varying the «j. 

Although linear term mixes the phonon modes, it is use- 
ful to define the non-renormalized energy of the phonon 
modes as 
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,,_£(,ma„9 (<, + ») + £!(<< + ») + 



where we have somewhat arbitrarily assigned jr of the 
linear and mass energies to each mode. 

By setting the mass to zero, we can confirm the division 
of energies implied by ([8|): i.e. the numerical minimisa- 
tion does in fact divide the energies equally between the 
masslcss quarks and the beads. 



VII. RESULTS 

In this section we show some results arising from this 
model. The intention is not to give a detailed compar- 
ison with experiment (which would require inclusion of 
colour Coulomb and magnetic interactions) but rather to 
show how the results change as we increase the number of 
beads and move away from the exactly soluble massless 
case. 

It is possible to get an intuitive feeling for the phonon 
states by using © to reconstruct the eigenmodes. We 
show the results for the 10-bead case for 3 combinations 
of quark masses; the case with two massless quarks (fig 
([3]) ) , with one massless and one heavy quark and 2 equal 
mass heavy quarks (fig (JSJ)). (fig 
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FIG. 3: Eigenmodes for Massless quarks 

Note that the heavy quarks stay relatively fixed, while 
the massless quarks move as much as the beads, as one 
would expect. The corresponding non-relativistic calcu- 
lation shows simply a straight line for the ground state 
solution for the case with two massive quarks: in the rel- 
ativistic case the beads necessarily carry a considerable 
amount of the energy. 




i i i i . 
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FIG. 4: Eigenmodes for one massive, one massless quark 
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FIG. 5: Eigenmodes for two massive quarks 
The renomalised energies of the bound state are given 

by 

i 

where ef (rij, ij) is given by (|10|) in the massless case. The 
only excited states we consider are the l\ — l,Zj>i = 
state, which corresponds to a conventional orbitally ex- 
cited meson, and the l\ = 0, \<x = 1, ^>2 = 0, which corre- 
sponds to the first hybrid state. If the quarks are mass- 
less, the ground state lies at E° = \/6hca ~ lGeV ( this 



10 



would be lowered in a realistic model by the short-range 
colour interactions). Fig © shows the energies of the 
li = 1, liyi = bound state as a function of quark mass 
for equal mass quarks and a varying number of beads. 
Note that the mass is dominated by the quark mass. 

hybrid, plol 




FIG. 6: Total energies for states consisting of two massive 
quarks and a varying number of beads. The lowest curve is 
the ground state (identical for all Nb by the renormalization 
process) and the upper curves show the energy for the h = 1 

It is more sensible to plot the excitation energy (i.e. 
the total energy with the quark mass subtracted) , and in 
Fig we show the ground state energy and the exci- 
tation energy of the 2 states for equal-mass quarks and 
varying number of beads. Fig ((8|) is the corresponding- 
plot for one heavy quark and one massless quark. The 
energies are roughly independent of the number of beads. 
The energies typically lie lower as the number of beads 
increases. The hybrid state lies consistently at twice the 
excitation energy of the first orbitally excited state 



VIII. CONCLUSIONS 



2500 




2000 Mq 4000 



FIG. 7: Excitation energies for states consisting of 2 massive 
quarks and a varying number of beads 
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FIG. 8: Excitation energies for states consisting of one mas- 
sive and one massless quark and a varying number of beads 



We have derived a relativistic model for mesons in a 
bead approximation to a flux-tube model. The model 
is analytically soluble for massless quarks, and can be 
renormalizcd to give a sensible spectrum. The massive 
quark case can be solved approximately by a rapidly con- 
verging algorithm. There are, however, several unattrac- 
tive features. 

1. The algorithm based on (fT9|) implies that the 
amount of energy carried by the quarks decreases 
as the number of beads increases. This is artificial, 



but it is not clear whether it has any real physical 
consequence. 

2. It is somewhat difficult to include short-range cor- 
rections: the colour Coulomb can be incorporated 
by a comparatively minor variation, as indicated, 
but the (crucial) hyperfine interactions probably 
require the extension of (JSJ) to include a Dirac equa- 
tion for the quarks. 

3. It is difficult to consider radial excitations in the 
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massive quark case. 

4. The choice of the renormalized ground state energy 
is somewhat arbitrary. 

These (and other) problems will be addressed in a later 
paper. 
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APPENDIX A: NOTATION 

There are several different energies in the calculation. 
They can be enumerated as follows: 

1. e i is the energy of the i th object (quark or bead) 

2. if (rii, U) is the non-renormalized energy of the i th 
phonon mode, excited to an (n i; ^) state, in the 
N-particle approximation. 

3. if (rii,li) is the renormalized energy of the i th 
phonon mode in the N-particle approximation. 

N 

4. E N (rii, k) = if (rii, U) is the non-renormalized 

i=l 

energy of the bound state. 

5. E° = E N (0, 0) is the ground state energy for the 
meson. 

N 

6. E N (rik,lk) — { n iih) is the renormalized 

i=l 

excited-state energy: i.e. the mass of the meson. 



APPENDIX B: GENERIC VARIATIONAL 
SOLUTION 

In all the calculations, the variational value of the en- 
ergy takes the following form: 

E 2 =J2 dU (a*) d) + KEd-! 09 fc ) + 

do {pj,ai) + adi (p k , fj,j,ati) + \ J2 3° 

' k 



with the following interpretations: 

1. (3k are the wave-function parameters.: i.e. -0 (ui) = 

2. The <Xi are the energy separation parameters: see 

en 



3. d k _ 2 (on) is the KE term: e.g for the 2-bead case 
with no Coulomb interaction, 

<± a = (§ + *+) 

r2„ 



d_ 2 (on) 



h z c z (3 
a B V2 



1-) 



4. c?„! (j3k) is Coulomb correction: note that it con- 
tains E since it is the 4th component of a 4-vector. 

5. do (iij, ai) is the mass term: 

do (Hj,ati) = 1 

so vanishes for massless quarks, 

6. d\ (cii, [ij, 0k) is the linear potential term: 



di (ai,Hj,Pk) = ( — + — 



for 0-bcad and 



d 1 (a i , H ,P i )=(^f^ + ^fW 

> — / JV „ N 1/2 



, N 



for the TVs-bead case. 
7. d\ (ai) arises from square of the linear terms: 

(§ + *+) 

OtB 



d i («*) = o + l - )id 2 (on) 



for 1-bead case 

The merit of writing it in this form is that it is com- 
pletely general and has an immediate solution for mass- 
less quarks with no Coulomb interaction: 



d*_ 2 M = Af , d\ (on) = K- t ,dU = di = d\ = o 

gives 

% 2 /Atv- 1 CT 



(A 



A, 2 



The general solution for the non-linear equation set can 
be written 



Pk = 

= Hk (Pi,Pk 



1/4 



(B2) 

This is the basis of the algorithm for rapid solution of 
the general case: 



1. Start with the above 



i0 



4 O') 

dt 2 (ai) 



1/4 
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2. Fix all Pi except for one, 0k- 

3. Find 3 successive values of Pk using eqn. (|B2|) 

4. Find the Aiken extrapolation to provide an up- 
dated [3k 



Pk 



P { k 2) Pk 0) -{Pl 1] ) 2 



5. Choose another Py and iterate until the values con- 
verge (typically 3 iterations for the full set). 

6. Use the final values of the Pi to find the energy E 
using IB 1 1 

Ths works very rapidly because 

Pjddx (aj.n^Pi) = 3 / /x t e 2 t fi N e 2 N \ 
dPk 2\a 1 fW + a N fWJ 

is a comparatively slowly varying function of the Pi. 

There is an extra complication when the Coulomb term 
is included because the energy needs to be calculated 
after each of the pi are estimated. Note that this format 
of the equation retains the phonon picture. 



APPENDIX C: STANDARD INTEGRALS 

For the one-bead case, we need to define a standard 
integral 



al ! i iC_ < ■ " - \u\ + ut\ k/2 u 2 + utdu + du_ 



r r 2Z + -8 2 u 2 2Z_ -R 2 U 2 

J J U + +p P+ U +m p P - U 



J 



e -y 2 {Pl cos 2 (0) +/ 3 2 sin 2 (0)) 

J cos 2+2Z + (0)sin 2+2i - (q>)y 5 + 2l ++ 2l -+ k dyd<j) 
r(3+z + +z_ + |) 

2 

2 + 2i + 2-\-2l_ 

cos (0)sin (0) 



(/3 2 cos 2 (0)+/3 2 sin 2 



since the y integral is trivial. Then the k = term is just 
the normalization: 

i _ O o (g 3 , _ r(f + * + )r(f + Z-) 



so that 



k (R R \ _ 9l + l_(P+,P-) 

~ A 2 

2 + 2l + 2 + 2l_ 

cos (4>) sin (0) 



ft +l _ (P+,P-) 



(01 cos 2 (0)+/3 2 sin 2 (4>)) 3+ ' + +l - +2 
This cannot be evaluated analytically for for k odd. How- 
ever, checks that can be run include: 



so approximately 
//,_ (/3+,/3_) 



'('+ + §) , ( L + §)' 



1/2 



/? 2 



APPENDIX D: N-BODY INTERACTIONS 

The trick in the previous section does not work for an 
iVs-bead solution with massive quarks, Instead we need 
to evaluate the expectation values of terms like 



1_ \ 



E 



ej e fc UjUk 



— /3 2 u 2 

We cannot evaluate / |ri| e * *Tlu 2 dui for k = 1 di- 
rectly. Write 



1/2 



JV n 2 

i— n+1 i— 1 1 



Ue-^u 2 du t 



(C is the normalisation) . Then we can expand, integrate 
and contract: 



2 u 2 



1/2 



^{A 1 +B 1 ) 1/2 = {A 2 + B 1 +e 2 u 2 ) 1/2 
or in general 



L4„ + S n _i + e„w 



Repeating this gives 



2,,2 I !/ 2 e -/3 2 « 2 ,2 



u 2 du n « (A„ + B n ) 1/2 



N 



1/2 



^ p 2 

=n-l ^ 



. All the cross terms vanish on integration. 

We can extend this: e.g. Coulomb term only exists 
between quarks, so write 



1 _ 1 

\fl— ?n] r 2 +r 2 — f 3 + ... — rjv| 

I J 



-1/2 



En 



jl ik 



giving the expectation value 

(Z N bcV /2 



E e H 
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